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ABSTRACT 

The large eccentricities of many giant extrasolar planets may represent the endpoint of gravitational 
scattering in initially more crowded systems. If so, the early evolution of the giant planets is likely to 
be more restrictive of terrestrial planet formation than would be inferred from the current, dynamically 
quiescent, configurations. Here, we study statistically the extent of the anti-correlation between giant 
planets and terrestrial planets expected in a scattering model. We use marginally stable systems 
of three giant planets, with a realistic range of planetary masses, as a simple model for the initial 
conditions prior to scattering, and show that after scattering the surviving planets reproduce well 
the known extrasolar planet eccentricities beyond a > 0.5 AU. By tracking the minimum periastron 
values of all planets during the evolution, we derive the distribution of orbital radii across which 
strong perturbations (from crossing orbits) are likely to affect low mass planet formation. We find that 
scattering affects inner planet formation at orbital separations less than 50% of the final periastron 
distance, (j/m, of the innermost massive planet in approximately 30% of the realizations, and can 
occasionally influence planet formation at orbital separations less than 20% of qfm- The domain 
of influence of the scattering massive planets increases as the mass differential between the massive 
planets decreases. Observational study of the correlation between massive and terrestrial extrasolar 
planets in the same system has the potential to constrain the origin of planetary eccentricity. 

Subject headings: solar system: formation — planets and satellites: formation — planetary systems: 
formation — celestial mechanics 



1. INTRODUCTION 

Studies of the population of small bodies - comets, as- 
teroids, and Kuiper belt objects - have provided pow- 
erful constraints on the early dynamical evolution of 
the massive planets in the Solar System. The concen- 
tration of objects in the outer Solar System in nar- 
row zones in resonance with Neptune, with significant 
eccentricities, provides evidence for an early expansion 
of Neptune's orbit (Malhotra 1995; Hahn & Malhotr^ 
ll999HMurrav-Clav fc Chianell2005|) . and has motivated 
a scenario in which the orbits of Saturn, Uranus 
and Neptune all expanded as a consequence of scat- 
tering the planetesim al debris left over after planet 
formation (Thommes. Duncan fc LevisonI 119991 120021: 
iTsiganis et al.l 20051. Migration of Jupiter - albeit for 
a much smaller distance - is also predicted within such 
models and may have left its own signature in the distri- 
bution of Hilda asteroids IjFranklin et al.ll200^ . 

There is no near or medium term prospect for ob- 
taining observations that would permit directly com- 
parable studies in extrasolar planetary systems. A 
wealth of circumstantial evidence, however, suggests that 
the early evolution of many extrasolar planetary sys- 
tems may have been more stochastic than in the So- 
lar System, and these more dramatic histories can be 
investigated by studying the relationship between ter- 
restrial planets (which are here being utilized purely 
as test particles of negligible mass) and giant plan- 
ets within the same system. The space-based tran- 
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sit mission Kepler IjBasri. Borucki fc Kochll2005D will of- 
fer the first definite opportunity to conduct such in- 
vestigations, though terrestrial planets may also be 
discovered via gi a nt planet transit t i ming techniques 
l|Agol et alJ I20M iHolman fc Murravl I200I or high- 
preci sion astrometry with the Space Interferometry Mis- 
sion (I Pord fc Tremaina 1200.111 ■ Gravitational microlens- 
ing also has potential for terrestrial planet discovery 
IjBennett fc Rhid l2002() . although the prospects for de- 
termining the overall architecture of systems in which 
planets are discovered are poor with this method. 

High eccentricities have become a common sig- 
nature of known extrasolar planets which reside 
outward of the tidal circularization limit. The ori- 
gin of thes e ecce ntricities remai ns unclear. Further, 
iGoldreichTL ithwic k fc Saril l|2004|) have suggested that 
the giant planets in the Solar System might have experi- 
enced large eccentricities before becoming circularized to 
their present values. One popular model attributes the 
eccentricity to gr avitational scatterin g among multiple 
massive planets llRasio fc FordI Il99(il:_ iLin fc Idal Il997t 
iLevison. Lissauer fc Duncan' 1998; 'Ford. Rasio fc Yul 
g003; Marzari fc Weidenschilling 2002). Scattering 
simulations have produced a quantitative match to the 
distribution of known e xtrasolar planet eccentricities 
IjFord. Rasio fc Yul |2003|) . and can also explain the 
current configuration of the upsilon Andromedae system 
|Ford. Lystad fc Rasio 2005). However, we cannot rule 
out the possibility that protoplanetary disk physics 
allows planetary eccentricity to be excited via the inter- 
action between a swqte massive planet and a gas disk 
JArtvmowicz e t a[\ 11991 ': 'Paoaloizou. Nelson fc MasseiJ 
,2001; Ogilvie fc Lubowl l2003: Goldrcich & Sari 2003). 

These two models for eccentricity growth have differ- 
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ent consequences for the formation of terrestrial plan- 
ets within the same system. If eccentricity arises from 
planet-disk interactions, then an eccentricity compara- 
ble to the final value is likely to be established soon after 
the formation of the giant planet, which must occur prior 
to th e dissipation of the gas dis k during the first few 
Myr IjHaisch. Lada fc Ladall2001() . This schedule prob- 
ably precedes the final stages of terr estrial p lanet for- 
mation, [chambers & Casscn (2002), ILevison & Agnor 
(|2003f) and IMardling fc LinI H2004(l have studied in de- 
tail the effect of the giant planet configuration, and its 
evolution, on the formation of terrestrial planets. If, con- 
versely, eccentricity is the endpoint of gravitational scat- 
tering among several massive planets, then the dynam- 
ics involved in the scattering process will usually involve 
multiple rearrangements of the planets prior to stabiliza- 
tion (usually by ejection) and establishment of a 'final' 
configuration. The extent of the ejection is dependent 
on the masses, radii, and initial orbital separations of 
the planets. During such rearrangements, which can take 
many Myr, giant planets can temporarily intrude (either 
physically, or via their resonances) into the terrestrial 
planet zone. The resulting constraints on where terres- 
trial planets can form and survive are bound to be more 
restrictive than those calculated assuming that the cur- 
rent observed (i.e. final) configuration of niassive planets 
has existed for all tim e (Jones. Slceo & Chambers 2001; 
i Menou fc TabachnikI 12003; .Asghari et aL.2004 .Ji et alt 

l2no5ir ^_^ 

In an earlier paper IjVeras fc Armitagel200{TD . we inves- 
tigated whether massive planet scattering - while pro- 
ducing final dynamically settled states that admit sta- 
ble terrestrial planets - might prevent terrestrial planets 
from surviving the transition to those states. Using test 
particles to follow explicitly the evolution of material (ei- 
ther formed terrestrial planets, or planetary embryos) in 
the terrestrial planet zone, we concluded from a limited 
number of realizations that there were circumstances in 
which the scattering process would impede subsequent 
planet formation at small radii. The dominant process 
for removing test particles was physical intrusion of a 
giant planet into the terrestrial planet region - i.e. the 
temporary existence of crossing orbits. 

In this paper, we quantify statistically the extent of the 
anti-correlation which massive planet scattering should 
introduce between the presence of massive and terrestrial 
planets in the sa me system. We assume, based on our 
previous results <|Veras fc Armitaggl200,'T|) . that terres- 
trial planet formation will be strongly suppressed exte- 
rior to the minimum periastron distance attained by any 
giant planet during the dynamical evolution, and calcu- 
late this minimum periastron distance from an ensemble 
of 500 N-body simulations. The appropriate initial con- 
ditions for scattering experiments are uncertain - here 
we assume that the initial states are triple planet point 
mass syste ms in which the planets have a realistic range 
of masses iTabachnik fc Tremain c""2002"). Although un- 
doubtedly oversimplified, these initial conditions lead to 
final states which match the observationally determined 
eccentricity distribution of exoplanets reasonably well. 

The outline of this paper is as follows. Section §2 de- 
scribes the motivation behind our initial conditions, and 
section §3 briefly describes the numerical codes used. We 
describe our results in §4, and discuss some of the impli- 



cations in S5. 



2. INITIAL CONDITIONS 



A large number of different initial configurations of 
massive planets are likely to evolve, via gravitational 
scattering, toward stable arrangements that match the 
observed eccentricity distribution of extrasolar planets. 
Massive planet formation theory does not provide any 
clear guidance as to the expected distribution of the num- 
ber of planets that would be formed in a single system, 
so it is reasonable to adopt any numbe r from N = 2 
{Ras io fc FordlllQQfiD to large N tPapaloiz OTi fc Tera^ieii] 
I2001D . That said, for simulations that ignore gas (as ours 
do), it is desirable for internal consistency to choose ini- 
tial conditions that yield instability on a timescale com- 
parable to, or longer than, th e dispersal timescale of the 
gaseous protoplanetary disk. iWolk fc Waited <)1996(1 es- 
timate this timescale observationally to be < 10^ yr, al- 
though it is not well constrained. Unfortunately many 
setups do in fact yield a significant fraction of rapidly un- 
stable initial conditions (for example Holman & Wisdoig 
|i993) find a logarithmic decay of the number of surviv- 
ing test particles in an outer Solar System simulation), 
though it is possible to at least avoid extremely short 
(and thus grossly inconsistent) instability timescales. 
Here, we model systems of three giant planets, with 
the motivation being that they represent the simplest 
systems which yield instability on the desired timescale 
without considerable fine-tuning of the initial planetary 
separations. 

2.1. Initial planetary separations 

iChambers. Wetherill fc Boss! l|1996(l and 

iMarzari fc Weidenschillina l|2002^ plot dynamical 
settling times for three-planet systems as functions of 
initial separations, while tito fc Tanikawa (2001) show 
similar plots for sets of 9 and 14 terrestrial-like planets 
which are perturbed by Jupiter and Saturn. The in- 
stability timescales for the three-planet systems exhibit 
an exponential dependence on initial separations, but 
were based on constant initial planetary separations 
and equal mass planets, simplifications not assumed 
in this work. Hence, although we used such diagrams 
as guidelines for establishing our initial conditions, we 
performed trial simulations of in order to determine the 
appropriate instability times for our runs. Based on 
those runs, we fixed our initial inner, middle, and outer 
planet semimajor axes at oi = 5 AU, 02 = 6 AU, and 
as = 12 AU respectively. The ratio of semimajor axes 
for the inner pair of planets, 02/01 = 1.2, lies just within 
the globally chaotic limit established by Wisdom (1980) 
for two-planet systems. Typically, these semimajor 
axis choices produce instability that causes the ejection 
of at least one planet on timescales exceeding several 
orbits but less than several million years. The system 
dynamically settles to a state that at least qualitatively 
resembles known multiplanet exosystems - often one 
featuring two well-separated eccentric giant planets. 

Since we ignore physical collisions between giant plan- 
ets, our simulations are freely scalable - only the ratios, 
and not the magnitudes, of 01,02,03 have physical sig- 
nificance. Nevertheless, our choices are plausible, since 
giant planet formation close (within a few AU) to a star 
is unlikely, while the predicted timescale for planet for- 
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mation via core accretion IjPollack et alJlTOOq) does not 
vary dramatically for a sa 5 — 10 AU. 

2.2. Masses 

Simulations of gravitational scattering within two- 
planet systems show that the re sults are dependent 
on the assumed mass d istribution IjRasio fc Fordl[l99a: 
iFord. Rasio fc Yull2003|) . Here, we randomly assigned 
masses to each giant planet in the range from O.SMjup to 
3Mjup, wei ghed by the 1/M power- l aw d istribution de- 
termined bv lTabachnik fc Tremaing 1)20021 . This range, 
the lower bound of which represents a Saturn mass, cov- 
ers one order of magnitude, and accounts for the major- 
ity of exoplanet minimum masses so far discovered. Our 
mass distribution is weighted toward subjovian masses 
in order to mirror the actual distribution of exoplanet 
masses. We note that we do not account for possible 
systematic differences between the mass distribution at 
formation, and that determined observationally following 
(within this model) gravitational scattering, which typi- 
cally ejects lower mass planets more frequently than high 
mass planets. Nor do we consider correlations between 
semimajor axis and mass at the epoch of formation. 

2.3. Eccentricities and inclinations 

Having determined the initial separations and masses 
of the giant planets, we randomized each planet's an- 
gular variables for each run. Each planet initially re- 
sides on a near-circular orbit, with an eccentricity of 0.01, 
and a random and negligible, but nonzero, inclination, of 
less than 0.001°. Such near-zero eccen tricity and inclina- 
tion va lues are similar to those used bv lFord. Rasio fc Yul 
()2003j) in their two-planet simulations. We fix our initial 
values so as not to introduce additional variables which 
should be considered w hen analyzing the r esults . Re- 
gardless, the variance in iFord. Rasio fc Yul 1)2003)) 's ini- 
tial eccentricity and inclination values do not appear to 
drastically alter their results. 

3. NUMERICAL METHOD 

We ran the scattering experiments using two inde- 
pendent gravitational N-body codes. Our primary tool 
was the HNbody code (K. P. Ranch fc D. P. Hamil- 
ton 2006, in preparation). HNBody is a collection 
of integrators primarily designed for systems with one 
massive central object, and has b een used for both 
massive e xtras olar planet simulations l)Veras fc Armitagg 
l2004albl 1200,'tD and satellite in teractions around Nep- 
tune ijZhang fc Hamilton! 12005)) . As a check, we also 
ran experiments with a Hermite stellar dynamics code 
(JHut fc Makino 2003), which, although not as efficient 
as HNbody for this problem, is perfectly able to evolve 
systems for the relatively modest durations of our sim- 
ulations. Adjustable run-time parameters for HNbody 
include choice of integrator, accuracy parameter, and 
timestep, whereas only timestep is a required parame- 
ter for the Hermite code. Given the Bulirsch-Stoer in- 
tegrator's ability to model accurately close encounters 
between planets, we chose this integrator for all HNbody 
simulations. The user-inputted timestep then sets the 
maximum timestep over which the integrator may iterate 
while still maintaining a prc-spccified accuracy. We used 
a timestep of 0.025 yr, which corresponds to a robust 
1/200 of an orbit of the initially innermost planet (note 
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Fig. 1. — The cumulative distribution of the observed extrasolar 
planet eccentricities beyond a > 0.5 AU (dotted lines) and the 
simulated distribution from the HNbody (solid lines) and Hermite 
(dashed lines) codes for 5 Myr simulations (upper panel) and 10 
Myr simulations (lower panel). The dot-dashed lines represents an 
additional simulated HNbody distribution derived from a different 
set of initial EGP semimajor axes of 5, 10, and 12 AU. 



that during the evolution, a for the innermost planet can 
decrease substantially, so a large safety margin at the 
start of the evolution is necessary). The resulting rela- 
tive energy errors were '~ 10~^ — 10~^^ throughout the 
simulation. For the Hermite code, we chose a smaller 
timestep of 1/1000 of an orbit for the innermost planet, 
as the code is not specialized for planetary problems. 
Even with the reduction in timestep, the relative errors 
(~ 10~* — 10~^) are orders of magnitude larger than 
those from the HNbody code, but still acceptable. Data 
was outputted at every 250-500 yr depending on the sim- 
ulation. 

In order to maximize the number of relevant simula- 
tions run while maintaining acceptable standards of en- 
ergy conservation, we performed a suite of preliminary 
tests. As the chaotic nature of the unstable systems 
we analyzed produces different evolutionary tracks due 
to machine-precision roundoff error, we do not compare 
simulations run with both codes on an individual basis. 
Rather, we perform a statistical comparison. We are also 
interested in comparing the distribution of final plane- 
tary eccentricities with observations, and for this purpose 
we consider only exoplanets with semimajor axes > 0.5 
AU; closer planets could have been influenced by tidal 
effects, which were not modeled by our simulations. We 
also only compare planets in systems in which dynamical 
instability produced crossing orbits, leaving only one or 
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TABLE 1 

KS TEST NUMBERS FOR ECCENTRICITY 

DISTRIBUTIONS 



5 Myr runs 



10 Myr runs 



actual-HNbody#l 

actual-Hermite 

HNbody#l-Hermite 

HNbody#2-HNbody#l 

HNbody#2-Hermite 

HNbody#2-actual 



13.8% 
71.5% 
41.0% 
88.3% 
75.0% 
75.6% 



36.0% 
57.6% 
46.6% 
96.0% 
22.9% 
41.2% 



Note. — The KS test probability that different pairs of 
distributions, plotted in Figure 1, are actually drawn from 
the same parent distribution. The designation "HNbody #2" 
represents the additional simulation run with initial EGP 
semimajor axes of 5, 10, and 12 AU. All four distributions 
are broadly consistent. 



two planets in stable orbits. Over half of the systems 
modeled with our initial conditions produced such insta- 
bility. 

4. RESULTS 

For our main sample, we evolved 250 systems us- 
ing HNbody for 5 Myr. This duration is compa- 
rabl e (in terms of orbits) to several previous stud- 
ies llFord. Havlickova fc Rasidl200H lAdams fc Lauehlinl 
120031 iVeras fc Armitaed l2004aD . Of these runs, 243 
achieved crossing orbits, and 175 exhibited at least one 
planet achieving a hyperbolic orbit. We also ran 250 Her- 
mite 5- Myr simulations, 244 of which achieved crossing 
orbits, and 170 of which exhibited at least one planet 
achieving hyperbolic orbits. No planet which remained 
stable throughout our simulations evolved its semima- 
jor axis to within a region where tidal effects would be- 
come important. Therefore, we did not need to incor- 
porate tidal corrections into the statistical comparisons. 
Although with our initial conditions significant dynam- 
ical evolution frequently occurs within 5 Myr, it is by 
no means obvious that the ultimate dynamical state is 
attained within this period. Hence, we carried out a 
smaller set of 50 10-Myr simulations with each integra- 
tor in order to check whether this alters the qualitative 
or quantitative results; we present results from both sets 
of simulations concurrently. The planets in all such simu- 
lations achieved crossing orbits. The number of systems 
which achieved hyperbolic orbits for the HNbody and 
Hermite simulations were 33 and 43 respectively. We 
only consider systems in which one or more planet has 
been ejected. 

Figure 1 displays the 91 known extrasolar planet ec- 
centricities as of June 1st, 2005 which satisfy a > 0.5 
(dotted line), along with the final eccentricity distribu- 
tions of the surviving planets from the 5 Myr (upper 
panel) and 10 Myr (lower panel) HNbody (solid line) 
and Hermite (dashed line) simulations. The HNbody 
distributions from a different set of initial EGP condi- 
tions, where oi = 5 AU, a-z = 10 AU, and 03 — 12 
AU, are superimposed on both panels as dash-dot lines. 
Respectable agreement is obtained, both between the re- 
sults from the two numerical codes, and between the sim- 
ulations and the observed distribution of eccentricities of 
extrasolar planets. Quantitatively, a KS test (Table 1) 
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Fig. 2. — Evolution of the semimajor axis (solid lines), periastron 
(dotted lines), and apoastron (dashed lines) of the three planets for 
a single run. In this case, the planet with a2 = 6 AU at 4 = (mid- 
dle panel) was eliminated from the system. The eliminated EGP 
performs a nearly complete radial sweep of the system between 1.5 
and 2 Myr. 



verifies that the simulated and observed eccentricities are 
consistent with their being drawn from the same parent 
distribution. This agreement may be somewhat fortu- 
itous, especially at the high e end of the distribution, 
where the observed population may have been un deres- 
timated due to selection effects jlCummind I2004H . Less 
than 2% of the planets simulated in this work attain ec- 
centricities higher than 0.8, the upper eccentricity limit 
predicted by Ford, Rasio & Yu ( 2003) for the two-planet 
scattering model in which both planets initially reside 
on coplanar, circular orbits. The occasional higher eccen- 
tricity values found in this study may be attributed to the 
four-body nature of the system and the small (= 0.01) 
nonzero initial EGP eccentricity values. 

The agreement demonstrated in Figure 1 suggests that 
our simulations yield outcomes that are consistent with 
the dynamical properties of actual systems. Of course 
this agreement is in no way unique, and other scattering 
scenarios are known that produce similar results. For 
example, we obta i n rat her similar results to those of 
iFord. Rasio fc YiJ l)2fl03t) . despite substantially different 
initial conditions. A range of masses among the scatter- 
ing planets does appear to be important. 

Within our model, the effective sphere of influence of 
the giant planets on subsequent terrestrial planet for- 
mation is dictated by the radial excursions of the giant 
planets during their dynamical evolution. Since dynam- 
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Fig. 3. — The distribution of the ratio Imllfm after 5 Myr, where 
gm is the minimum periastron distance achieved by any planet 
throughout the simulation, and qf^-n is the minimum periastron 
among the surviving planets at the final state of the simulation. 
HNbody simulations are represented by solid lines and Hermite 
simulations are represented by dashed lines. 
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Fig. 5. — The distribution of the ratio q-m/qfm after 10 Myr, 
where qm is the minimum periastron distance achieved by any 
planet throughout the simulation, and qfm is the minimum pe- 
riastron among the surviving planets at the final state of the sim- 
ulation. HNbody simulations are represented by solid lines and 
Hermite simulations are represented by dashed lines. 
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Fig. 4. — The cumulative distribution of the ratio qm/qfm after 
5 Myr, where qm is the minimum periastron distance achieved by 
any planet throughout the simulation, and qjm is the minimum 
periastron among the surviving planets at the final state of the 
simulation. HNbody simulations are represented by solid lines and 
Hermite simulations are represented by dashed lines. 



ical instability often causes crossing orbits, the most re- 
strictive conditions (i.e. the greatest inward excursions) 
can be set by any of the initial three planets. One exam- 
ple is depicted in Figure 2, which shows the evolution of 
the semimajor axis, pericenter and apocenter of the three 
planets for the first 2 Myr of a 5 Myr run. In this case, 
the final state of the system leaves two planets in moder- 
ately eccentric, well-separated periodic orbits. However, 
in the evolution up to this state, the initially innermost 
and outermost planets switched positions, while the third 
planet swept radially through the entire system within 
2 Myr, and harbored a pericenter which remained close 
to the parent star for nearly 0.5 Myr. This behavior is 
effectively hidden by the benign final state of the system. 
In order to quantify how common such evolutionary 
histories are, we consider the ensemble of simulations 
discussed above. We denote by g„i the minimum peri- 
astron distance achieved by any planet throughout the 
simulation, and by qfm the minimum periastron among 



the surviving planets at the termination of the simula- 
tion. Figure 3 shows the distribution of the ratio of these 
distances, qm/qfm, for the 5 Myr HNbody (solid line) 
and Hermite (dashed line) simulations respectively. Fig- 
ure 4 is the correspondin g set of cumulative histo grams. 
We find, as previously IjVeras fc Armitaggl2005|) . that 
the periastron distance for the innermost planet in the 
final, relaxed state, is often a poor proxy for the extent 
of radial sweeping that occurs during the earlier dynam- 
ical evolution. Although qm/qfm does increase toward 
Qm/qfm = Ij the distribution is broad and extends down 
to quite small values qm/qfm ^ 0.2. From the cumula- 
tive histograms, we find that qm/qfm ^ 0.2 in approx- 
imately 10% of realizations, while qm/qfm ^ 0.5 in ap- 
proximately 30% of cases. Within a scattering scenario, 
we would therefore expect to observe a significant deficit 
of lower mass planets in orbits with semimajor axes less 
than half the periastron distance of the innermost mas- 
sive planet. Many orbits in this exclusion zone would be 
stable today, but would have been unfavorable for planet 
formation due to the early massive planet dynamics. No 
significant difference between results computed with the 
two codes is seen. Broadly similar results are also ob- 
tained from the simulations run on for 10 Myr - plotted 
in Figure 5 - albeit with poorer statistics. 

No simulations so far described include terrestrial plan- 
ets. As a further demonstration of our results and as a 
check on the approximate criterion for strong interac- 
tions adopted earlier, we performed 7 additional simu- 
lations including 50 terrestrial planets over 5 Myr with 
HNbody. The planets were treated as "Light Weight Par- 
ticles" (LWPs), which are test particles that do not per- 
turb EGPs. The LWPs were initially placed uniformly 
in annuli ranging from 0.5 — 1.5 AU (2 simulations) or 
1.0 — 3.0 AU (5 simulations) on circular, nearly copla- 
nar orbits. Figure 6 provides a snapshot of the final 
states of all 7 simulations, with dots representing surviv- 
ing LWPs, solid lines representing g^, and dotted lines 
enclosing the initial LWP annulus. The plots are ar- 
ranged according to decreasing g,„ in order to illustrate 
just how sensitively terrestrial planet survival depends 
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distance achieved by any EGP throughout the simulation. Two 
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upon EGP excursion into the inner system. These runs 
show that although oceasionaUy low mass bodies may 
survive with semi-major axes outside the minimum pe- 
riastron separation, typically there's significant clearing 
even within the minimum periastron distance. Hence, 
the approximate criteria we've adopted is on the conser- 
vative side. The uppermost panel indicates that the ini- 
tially low eccentricities of LWPs will remain near-circular 
if the EGPs fail to intrude upon a region within twice the 
radial distance of the terrestrial material; the seemingly 
anomalous vertical band of LWPs at « 1.3 AU may be 
attributed to strong resonant interactions with the ap- 
proaching EGPs. 

The random planet mass distribution used in the sim- 
ulations raises an additional question - which combina- 
tions of planet masses promote the greatest extent of 
radial sweeping? Under some circumstances (for exam- 
ple when two planets remain on stable bound orbits at 
the end of the evolution), clues as to the initial masses 
of the scattering planets may survive the scattering pro- 
cess, so any strong correlation between the masses and 
the outcome is potentially observable. Figure 7 shows 
how the mean q^ scales with the range of planet masses 
in the randomly chosen initial conditions, quantified via 
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Fig. 7. — Variation of the mean qm values with the variance 
of the initial planet masses, for 5 Myr (solid line) and 10 Myr 
(dashed line) simulations using the HNbody (triangles) and Her- 
mite (squares) codes. Plotted points refer to binned values of the 
logged mass standard deviation. Masses were randomly sampled 
from a 1/M distribution. 



the standard deviation of the log of the planet masses. 
In the figure, solid lines denote the 5 Myr simulations, 
dashed lines indicate 10 Myr simulations, and triangles 
and squares refer, respectively, to the HNbody and Her- 
mite code simulations. We find a roughly linear trend - 
as the scatter in the mass increases, so does the mean 
value of q^n ■ Systems in which the planets are of roughly 
similar masses are therefore most detrimental to the sur- 
vival of interior terrestrial planet material. Given our 
assumed 1/M distribution of masses, we expect EGP 
standard mass deviations of over IMjup to be relatively 
rare. 

Within our assumed scattering scenario for the origin 
of extrasolar planet eccentricities, the extent of the pen- 
etration of giant planets into the terrestrial planet zone 
is likely to be the dominant factor determining whether 
the evolution hinders subsequent terrestrial planet for- 
mation. However, the amount of time spent with giant 
planets at small periastron distances may also be signifi- 
cant. We find a wide gap between the timescale on which 
crossing orbits first develop, and the timescale on which 
a planet is first forced onto a hyperbolic orbit. Crossing 
orbit times for our simulations are on the order of 10"^ 
yr - short enough that the effects of a residual decaying 
gas disk would obviously be important and modify the 
result. By contrast, the timescale for ejections is of the 
order of a Myr. Figure 8 is a histogram of the initial 
time on which a planet achieves a hyperbolic orbit in the 
5 Myr simulations for the HNbody (solid line) and Her- 
mite (dashed line) simulations. Because the timescales 
on which planetary orbits merely cross versus become 
hyperbolic vary by several orders of magnitude, Saturn 
or Jovian-mass planets may be dynamically excited for 
hundreds of thousands of years before being ejected from 
the system, suffering a collision, or dynamically settling 
into a quasi-stable configuration such as that in Figure 2. 
Intrusions into the terrestrial planet zone are therefore 
typically repeated and / or extended, rather than being 
one-time events. 

5. DISCUSSION 
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Fig. 8. — Histogram of times at which a planet first traveled on 
a hyperbolic orbit. The planets were evolved for 5 Myr. HNbody 
simulations are represented by solid lines and Hermite simulations 
are represented by dashed lines. 



Forthcoming missions that promise to detect terres- 
trial planets within the habitable zone - such as Ke- 
pler - will also open the possibility of studying sta- 
tistically the influence of massive planets on terres- 
trial planet formation (at least in the limited cases 
in which radial velocity surveys prove sensitive enough 
to constrain the giant planet populations around typi- 
cally rather faint stars). That the existence of Jupiter 
and Saturn affects the formation of Earth and the rest 
of the inner Solar System, particularly in the vicin- 
ity of the asteroid belt, is well-established llWetherilll 
r|9^; "Lecar & Frankhn"199r; 'Chambers fc Casscnl l200'a 
Eevison & Agnor 2003: Mardling & Lin 2004). The dis- 
tribution of extrasolar planetary masse s is securely pre- 
dicted to show a similar bimodality ijida fc LinI 120041 
e.g.), due to the rapid timescale on whic h runaway ac- 
cretio n of giant planet envelopes occurs IjPollack et al.l 
Il996(l . so a coupling between the observed properties of 
the giant and terrestrial planets in the same system is ex- 
pected. Since giant extrasolar planets have substantially 
different orbital properties (and, perhaps, evolutionary 
histories) than those planets in the Solar System, we 
have argued here that the coupling between giant and 
terrestrial planet formation might be atypically weak in 
the Solar System, and more dramatic in extrasolar plan- 
etary systems. 

The orbital properties of giant extrasolar planets 
immediately suggest two new mechanisms via which 
the evolution of giant planets might impede terres- 
trial planet formation. First, the orbital migration 
that is believed to b e necessary to explain the ori- 

tof Hot Jupiters l)Lin. Bodenheimer fc RichardsonI 
7), and possibly the Late Heavy Bombardment 
)mes et al. I '20051. means that, at le ast in some sys- 
tems and possibly in the typical system IjArmitaee et al.l 



l2002HTrilling. Lunine fc Bendl2002|) . giant planets move 
through the entire terrestrial planet zone during the early 
evolution of the disk. This may reduce the probability of 
subsequent terrestrial planet formation (Armitagc 2003), 
though the extent of any suppression will depend on the 
speed of migration (which alters the degree of shepherd- 
ing via resonant capture) and on the surface density of 
planetesimals in the i nner reg ions following scattering by 
the migrating planet llWard fc Hahn 1995: Tanaka fc Id^ 



1999; MandeU fc Sigurdsson 2003; Edgar fc Artvmowicz 
2004; Ravmon d. Quinn fc L unine 2005; Fogg fc Nelso' 
2(}05). EGPs which accrete onto the star may elongate 
the dynamical relaxation time of the planetesimal density 
distribution, and may occur repeatedly, further inhibit- 
ing terrestrial planet formation. Second, the large mean 
eccentricities of known giant extrasolar planets may sig- 
nal that the typical planetary system is initially dynam- 
ically crowded and unstable to eventual gravitational 
scattering and planetary ejection. If this is the origin of 
eccentricity, then transient incursions by the giant plan- 
ets into the terrestrial planet during the scattering pro- 
cess can be particularly detrimental to terrestrial planet 
formation, since the scattering is likely to occur at a rel- 
atively late epoch when the gas disk (which would tend 
to recircularize scattered planetesimals) has been dissi- 
pated. Studies of the dynamical stability of hypothetical 
terrestrial planets in known extrasolar planetary systems, 
which are normally based on presently observed pericen- 
ter values, may then hide the extent to which planets 
have swept out terrestrial material earlier in the history 
of that system ( Vcras fc Armitagc 2005) . 

In this paper, we have attempted to quantify - within a 
specific scattering scenario - the extent (in radius) of the 
anti-correlation this mechanism would introduce between 
the presence of giant and terrestrial planets in the same 
system. We find that during evolutionary sequences that 
end in plausible-looking planetary systems, one massive 
planet sweeps inward to less than half of the final ('sta- 
ble') periastron distance approximately 30% of the time. 
Some effect persists (at the 10% level) down to orbital 
radii as small as 20% of the final periastron distance. 
If these simulations reflect the physical processes tak- 
ing place in the early history of real extrasolar planetary 
systems, we predict that terrestrial planets ought to be 
under-abundant in nominally dynamically stable orbits 
interior to eccentric massive planets. 
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